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Rice fragrance is an important characteristic for Southeast Asian consumers, and fragrant landraces from 
Japan were first recorded in the 17th century. Principal component analysis clearly showed that Japanese fra- 
grant landraces were genetically different from non- Japanese fragrant landraces. Japanese fragrant landraces 
were composed of six clades, none of which carried the most common fragrance mutation, an 8-bp deletion in 
exon 7 of Badh2. Fragrant landraces comprised two major groups carrying different Badh2 mutations. One 
group carried a known SNP at exonl3 and the other a SNP at the exonl-intronl junction as splicing donor site. 
The latter was considered to be a potential splicing mutant group as a novel allele at Badh2. Heterozygosity 
{He) scores in the two fragrant groups were not significantly different from non-fragrant landraces and modern 
cultivars. However, lower He scores were found around the Badhl locus in the two groups. The potential 
splicing mutant group showed a more extended haplotype than the E13 SNP group. A likely causal factor re- 
sponsible for loss of function is a novel splicing mutation allele that may have been generated quite recently. 
The fragrance allele has dispersed as a result of out-crossing under local environmental conditions. 

Key Words: fragrance, landrace, 2-acetyl pyrroline, Badh2, diversity, selective sweep. 



Introduction 

Rice consumers in Southeast and South Asian countries 
have a preference for fragrant rice, which has been selected 
from fragrant landraces that were initially improved by se- 
lection of pure lines from natural variants. Fragrant land- 
races have been considered a valuable resource in comparison 
with modern cultivars, and many genetic analyses of these 
landraces have been reported (Ahn et al. 1992, Amarawathi 
et al. 2008, Bradbury et al. 2005, Bourgis et al. 2008, Chen 
et al. 2006, Lanceras et al. 2000, Lorieux et al. 1996, Shi et 
al. 2008). Fragrance has been reported to be regulated by a 
single or multiple genes (Berner and Hoff 1986, Huang et 
al., 2008, Sood and Siddiq 1978) and although a major 
causal gene, Badh2, has been located on chromosome 8, 
other loci have also been reported to affect rice fragrance 
(Amarawathi et al. 2008, Lorieux et al. 1996). Loss of func- 
tion of Badh2 results in accumulation of 2-actyl-l -pyrroline 
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(2AP) (Bradbury et al. 2005, Vanavichit and Yoshihashi 
2012), the major component of rice fragrance (Magal984, 
Paule and Powers, 1989, Yajima et al. 1978). The most 
widely distributed polymorphism in Badh2 that causes 2AP 
accumulation is an 8-bp deletion in exon 7. Up to now, ten 
alleles including the major mutation have been reported 
(Kovach et al. 2009, Shi et al. 2008), all of which originate 
from japonica and have introgressed into other genetic 
backgrounds. 

'Kabashiko' is the generic name for fragrant rice in 
Japan, and various fragrant Japanese landraces primarily 
distributed in mountainous areas, mentioned in 17th century 
historical records, Aizu-nousho, a document published in 
1684 providing a record of agricultural information (Arashi 
1975). By the end of the 19th century, application of strin- 
gent quality control had excluded red rice, slender grains, 
low-quality grains and so on. Most of the fragrant landraces 
also disappeared around that time. Why, then, was fragrance 
considered undesirable? Fragrant landraces may be exluded 
probably because of other undesirable traits such as a lower 
panicle yield (Arashi 1975). It is also true that the full un- 
diluted smell of fragrant grains would probably have been 
unpopular in Japan. In general, consumers added a small 
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portion of fragrant grains to aged grains to add a fresh 
flavor. Only restricted usage was applied in Japan. How- 
ever, fragrant landraces were dispersed into Japan gradually 
during 18th to 19th centuries. These fragrant landraces 
show geographical divergence in their traits (Itani 2002a, 
Miyakawa and Nakamura 1984). Based on morphological 
and physiological traits, they can be roughly classified into 
two major groups: eastern and western. The eastern group 
distributed in Hokuriku, Kanto, and Tohoku, was character- 
ized by early heading, short panicles, low yield, a fine culm 
and a propensity to lodge, relative to the western group dis- 
tributed in Kinki, Shikoku, and Kyushu. Like red rice, they 
soon disappeared from the general market and the fields 
used for growing them were resown with modern cultivars 
over the last 100 years (Arashi 1975, Ishikawa et al. 2002a, 
2002b, Itani 2002a, Tabuchi et al, 2007). Although land- 
races including fragrant landraces were displaced by modern 
cultivars they may retain a number of potentially valuable 
traits that have been lost in modern cultivars (Arashi 1975). 

In the present study, we analyzed genetic diversity 
among Japanese fragrant landraces in comparison with 
Japanese non-fragrant landraces and modern cultivars. Ge- 
netic distance allowed us to specify particular groups of 
landraces that had adapted to local environmental condi- 
tions and carried a novel mutation related to fragrance. This 
novel mutation implied us different amount of fragrance 
based on the gene structure. Thus, it will allow us to design 
fragrance corresponding to needs of cunsumers' preference. 

Materials and Methods 

Plant materials 

Eighty-two cultivars, including landraces and modern 
varieties, were examined (Table 1). These landraces were 
collected at Hiroshima prefectural Hiroshima University 
and Hirosaki University. Most of those landraces were char- 
acterized for morphological and physiological traits (Itani 
2002a, 2002b). Of these, 17 were non- Japanese fragrant 
landraces, 25 were Japanese non-fragrant cultivars includ- 
ing elite cultivars such as 'Koshihikari', and 40 were 
Japanese fragrant landraces included five non-scented 
ones. Overseas category included Meragome originated 
in Kagoshima which is an indica cultivars known as 
Daitoumai-group introduced into Japan in the middle age 
(Ishikawa 2002a, Ogawa and Itani 2008). Thus, this land- 
race was categorized as a non-Japanese landrace. 
'Masshigura' and 'Koihonoka' were used to confirm allelic 
relations among various fragrant landraces. Both of them 
are improved varieties in Aomori prefecture. The former 
one is non-fragrant leading variety. The latter one is a fra- 
grant rice inherited 8bp deletion in Badh2 gene from 'Sari 
Queen'. 'Sari Queen' was generated from the cross between 
'Basmati370' and 'Nipponbare'. Then, the variety was 
crossed with other Japanese varieties to obtain 'Koihonoka' 
(Ando et al. 2004, Mikami et al. 2007). 



Fragrance testing by treatment with KOH 

The fragrance (2AP) status of cultivars was confirmed by 
soaking small pieces of leaves in 2% KOH solution heated 
to 60°C for 10 minutes. Sensory analysis was performed by 
multiple persons at multiple times to confirm phenotypes. 
The varieties were also classified as tropical or temperate 
japonica (Tr or Tm) using the method described by Sato 
(1991). 

Molecular markers 

Thirty SSR loci dispersed across the genome (Supple- 
mental Table 1) were amplified with Ampliqon-Taq 
(Ampliqon Co., Denmark). PCR conditions were 94°C for 
3 min, 30 cycles of 94°C for 10 s, 55°C for 30 s, and 72°C 
for 30 s, followed by 72°C for 5 min for post-heating. Am- 
plified products were genotyped using 6% denatured se- 
quencing gels with a general silver staining method based 
on the SILVER SEQUENCE™ DNA Sequencing System 
(Promega Co., USA). Additional SSRs were analyzed in 
Japanese landraces and elite cultivars to screen for a selec- 
tive sweep around Badh2 (Supplemental Table 2). INDEL 
markers were developed for Badh2 from representatives of 
fragrant landraces (Supplemental Table 3). The PCR proto- 
col was the same as that for SSRs. A CAPS marker includ- 
ing the splice donor site mutation at the exonl-intron 1 
junction was developed to confirm whether the SNP was 
cleaved by the restriction enzyme Kpnl (New England Bio- 
Labs, Japan, Inc.) using PCR fragments amplified with the 
5'-UTR forward primer and 1-986 reverse primer shown in 
Supplemental Table 3. 

Badh2 sequences 

Parts of Badh2 were amplified with primer pairs listed in 
Supplemental Table 4. Amplified fragments were purified 
with a QIAquick PCR purification kit (QIAGEN Co., 
Japan). These fragments were then prepared for sequencing 
using a BigDye Terminator v3.1 Cycle Sequencing kit 
(Applied Biosystems Co., Japan) to obtain sequences in- 
cluding the ABB 10 sequence. SNP mutation found at the 
exon-intron junction has been registered (AB821274). 

Tr-japonica, Tm-japonica classification 

Hull size, mesocotyl length, and KOH endosperm degra- 
dation were tested using a protocol described by Sato 
(1991) (Supplemental Table 5). Hull size was measured us- 
ing Nikon Digital Sight (NIKON, Japan) and expressed as 
the average for ten hulls. Mesocotyl length (mm) was mea- 
sured after two weeks of germination in total darkness. Ten 
young seedlings were subjected to 1st internode measure- 
ment and the values averaged. KOH degradation was mea- 
sured as the loss of weight after soaking in 1.5% KOH solu- 
tion, ten grains being thus treated. Z2 score was determined 
according to Sato (1991): 

Z2 score = R + 0.093M - 0.044E - 2.01 

where R: ratio of hull length to width, M: mesocotyl length 
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Table 1. List of materials characterized genotypes and phenotypic traits of landraces (LR) and modern varieties (MV) 



Overseas fragrant landraces 



Japanese fragrant landraces 



Japanese non-fragrant landraces/modern varieties 
(MV) 









Fra- 










Fra- 










Fra- 




Variety name 


Origin 


Criteria 


grance 


Clade* 


Variety name 


Origin 


Criteria 


grance 


Clade 


Variety name 


Origin 


Criteria 


grance 


Clade 


Carolina 


USA 


LR 


+ 


4 


Nioi-wasel 


Iwate 


LR 


- 


1 


Aichiasahi 


Aichi 


MV 


- 


6 


Kara lath 


India 


Lis. 




i 


Kaon-wase 


Iwate 


i \? 

LK 




1 
1 


Aikoku 


Shizuoka 


LK 




0 


Ikkikoku 




LR 


-l_ 


4 


Jakou-ine 3 


Yamag at a 


LR 




[ 


Akage ^ 


Hokkaido 


LR 






Meragome 


Kagoshima 


LK 




-i 

4 


Kaon-ine 5 




i 

LK 




1 


Akamai 3 


Unknown 


LK 




o 


Mu sk akudant 1 




LR 


_l_ 


4 


Kaon-wase 2 


riokJtaiQo 


LR 






Akihikari 


Aoniori 


MV 




(•-, 


S u rj amkh 1 


India 


LR 




4 


Kaon-wase 3 


Hokkaido 


LR 




2 


Akitakomachi 


Akita 


MV 




(•-, 


Choukokou 


China 


i d 




4 


Kaon-ine 4 




1 [.> 
LK 




z 


Asahi 


Kyoto 


LK 




0 


Azucena 


Philippines 


Lis. 




c 

j 


Hokkaido- 
40nichi-wase 


Iwate 


i \? 

LK 






Fujisaka 5 


Aomon 


1V1 V 




0 


rsasmati j /u 


India, 
Pakistan 


LK 




j 


Jakou (Niigata) 


Niigata 


1 D 

LK 


+ 


i 
1 


Ginbozu 


Toyama 


\ AX 7 
1V1 V 




0 


Bimcol 


Philippines 


1 D 

LK 


+ 


c 

J 


Nioi-wase2 


Akita 


LK 


+ 


i 
1 


Gohyaku- 
mangoku 


Niigata 


LK 




0 


... 

Beijmg-koutou 


China 


LR 


+ 


5 


Iwaka 


Miyagi 


LR 


+ 


3 


Hyogoomachi 


Hyogo 


MV 




6 


DawDam 


Thai 


LR 


+ 


5 


Jakou-ine 2 


Tohoku 


LR 


+ 


3 


Jukkoku 


Unknown 


LR 






Delia 


USA 


LR 


+ 


5 


Jakou (Ibaraki) 


Ibaraki 


LR 


+ 


3 


Jyoshu 


Yamagata 


LR 




6 


KliaoPickLaos 


Laos 


LR 


+ 


5 


Jakou-wase 


Niigta 


LR 


+ 


3 


Kameji 


onimane 


LR 




6 


KU701 


Thai 


LR 


+ 


5 


Kaori-ine 1 


Ibaraki 


LR 


+ 


3 


Kamenoo 


Yamagata 


LR 


- 


6 


Mangareez 


Philippines 


LR 


+ 


5 


Kaori-ine 2 


Kanagawa 


LR 


+ 


3 


Koshihikari 


Fukui 


MV 




6 


Sabarmati 


India 


i d 

LK 


+ 


c 

J 


Nioi-akako 


Miyagi 


LK 




i 

j 


Nipponbare 
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1V1 V 




/: 
0 
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Miyagi 


LK 
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1 I? 
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o 
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l D 

LK 






t? ; Vim i n 

K1KUU 1 JZ 
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1V1 V 
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koutou 






























Kum amoto 


Kumamoto 


1 Q 

LK 




1 


. ... . 
Sasanishiki 


Miyagi 


1V1 V 




0 












Kabashiko A 


Miyazaki 


LK 




1 


Senichi 


Tochigi 


LK 




o 












Kabashiko B 


Kuma mot o 


LR 
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J 


Shinnki 


Hyogo 


LR 




5 












JvdDdSniKO 
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LR 




7 


Takenan 


ivf 16 1 ' 


LR 
















Kaon-ine 3 


Nara 


LR 




7 


Tsushima 


JN agasaki 


LR 
















ICeitoku 


K.um amoto 


Lis. 




7 

1 


Tsushima- 
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Kumamoto 
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Kumamoto 
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7 






















Nioi-yoshi 1 


Kochi 


LR 


— 
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Nioi-yoshi 2 


Kochi 


LR 


— 


7 






















Seigoku- 


Nagasaki 


LR 


— 


7 






















kawachi 






























Shiro-wase 


Wakayama 


LR 


- 


7 






















Jakou-ine 1 


Yamagata 


LR 


+ 


8 






















Jakou-ine 4 


Yamagata 


LR 


+ 


8 






















Jakou-mochi 


Touhoku 


LR 


+ 


8 






















Kamari 


Ibaraki 


LR 


+ 


8 






















Nioi-mochi 1 


Fukuhsima 


LR 


+ 


8 






















Nioi-mochi 2 


Fukushima 


LR 


+ 


8 






















Hieri 


Kouchi 


LR 


+ 


9 






















Kaisen 


Kouchi 


LR 


+ 


9 






















Wase-hieri 


Kouchi 


LR 


+ 


9 













! Clade was based on genetic distance shown in Fig. 2. 

! * Landrace introduced from overseas was classified into oversea origin. 



(mm), E: reduction of endosperm weight after soaking in 
1.5% KOH. Cultivars showing a Z2 score of <0 were identi- 
fied as Tm-japonica (Tm), and those showing a score of >0 
as Tx-japonica (Tr). 

Data analysis 

The GenAlEx 6.2 software package (Peakall and Smouse 
2006) was used to evaluate genetic variation in loci and 
populations, including the number of different alleles per 
locus (Na), Fst, observed heterozygosity (Ho), and expected 
heterozygosity (He). Dendrograms were constructed using 
the neighbor-joining method based on Nei's unbiased genet- 
ic distances using the Populations 1.2.30 beta2 program 
(http://bioinformatics.Org/~tryphon/populations/#ancre_ 



bibliographic). All dendrograms were drawn by TreeExplorer 
(Tamura et al. 2011). Bootstrap values were assessed with 
1000 replicates. 

Results 

Diversity in fragrant rice 

Confirmation of fragrance found five Japanese fragrant 
landraces, 'Nioi-wase 1', 'Kaori-wase', 'Jakou-ine 3', 
'Kaori-ine 5', and 'Narukosan-koutou', revealed that they 
were not fragrant, even though they shared similar names 
implying fragrance (Table 1). All other Japanese and non- 
Japanese fragrant landraces were fragrant. Including these 
landraces, genetic divergence among landraces and modern 
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Principal Coordinates 
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Fig. 1. Principal component analysis of fragrant rice by using 30 
SSR loci. Overseas fragrant landraces (squares), Japanese fragrant 
landraces (diamonds), and Japanese non-fragrant landraces/elite varie- 
ties (Triangles) were separated clearly. 



cultivars was genotyped using 30 SSRs. 

Principle component analysis (PCA) clearly identified 
three major groups: overseas (non- Japanese) fragrant land- 
races, Japanese fragrant landraces, and Japanese non-fragrant 
landraces/varieties (Fig. 1). The non- Japanese fragrant land- 
races were the most diverse, with He = 0.652 (Table 2). 
These were followed by Japanese landraces and cultivars in 
combination, He = 0.397, Japanese non-fragrant landraces 
and modern cultivars with He — 0.349, and Japanese fra- 
grant landraces with He = 0.336. There was no significant 
difference between the Japanese fragrant and non-fragrant 
groups. 

Genetic relationships among fragrant landraces 

A phylogenetic tree was constructed based on genetic 
distance (Fig. 2). Overseas cultivars were grouped into 
Clades 4 and 5. Landraces in Clades 4-2 and 5 carried an 
8-bp deletion, which is known to be the major causal muta- 
tion in Badh2, but not in Clades 4-1 or in other Japanese 
fragrant landraces. Japanese non-fragrant landraces and 



modern varieties were included in Clade 6. Japanese fra- 
grant landraces were classed into six groups (Clade 1, 2, 3, 
7, 8, and 9). Clade 1 comprised fragrant landraces giving no 
aroma in the sensory test, but which had names related to frag- 
rance and originated from Iwate and Yamagata in Tohoku 
district. Varietal names similar to them, such as 'Kaori- 
wase' and 'Nioi-wase', were also found in Clade 2, and 
were phenotyped as fragrant in the sensory test. Landraces 
in Clade 2 originated from Eastern Japan, extending from 
Niigata to Hokkaido. Clade 3 included one fragrant land- 
race, 'Narukosan-koutou', which gave no aroma in the 
sensory test. In addition, it was most distantly separated 
from the other members of the clade. Fragrant landraces in 
Clade 3 originated from Tohoku area to Kanto area (Ibaraki 
Prefecture). Clade 7 comprised landraces from Western 
Japan. Clade 8 included three glutinous and three non- 
glutinous landraces, all of which originated from northern 
Kanto up to the southern areas of Tohoku. Clade 9 com- 
prised three landraces that originated from the same prefec- 
ture, Kochi, in Western Japan. 

Novel fragrant allele at the Badh2 locus 

In order to confirm molecular variation in Badh2, the 
full-length Badh2 gene including the 1296-nt upstream re- 
gion was sequenced from representatives of each clade 
(Fig. 3). A C to T (C/T) substitution as a single nucleotide 
polymorphism (SNP) in exon 13, leading to substitution of 
a single amino acid, has already been reported in a single 
fragrant landrace, 'Shiratama', originating from Japan. 
'Hokkaido-40-nichi-wase', 'Iwaka', 'Kabashiko A', and 'Hieri' 
belonging to Clades 2, 3, 7, and 9, respectively, were se- 
quenced and showed the C/T substitution, which results in a 
single amino acid change (alanine to valine) (Table 3). 
However, 'Nioi-mochi 2' in Clade 8 did not carry the muta- 
tion. Exon 13 was sequenced in all Japanese landraces, and 
all fragrant landraces carried the substitution at exon 13 ex- 
cept for those in Clade 8. Another significant nucleotide 



Table 2. Heterozygosity among groups with different genetic background such as overseas fragrant group carrying 8bp deletion or no-deletion at 
exon 8, non-scented group, E13 SNP group carrying a SNP in exon 13, splicing donor site mutant (splic-mut) group, and non-fragrant group 







He (all loci) 






He around Badhl 






Group 


No. of - 
cultivars 


Mean ± SE 


544kup 
(RM8264) 


45kup 


12kup 


SSRat7 


1 06kdwn 
(RM1109) 


142kdwn 
(RM3459) 


Total 


82 


0.516 ± 0.043 


0.736 






0.510 


0.668 


0.775 


(Total — Non-scented) 


77 


0.523 ± 0.043 


0.741 






0.528 


0.676 


0.771 


Oversea 


17 


0.652 ± 0.028 


0.761 






0.734 


0.699 


0.747 


8 bp-del. 


7 


0.467 ± 0.044 


0.612 






0.490 


0.408 


0.694 


non 8 bp-del. 


10 


0.523 ± 0.042 


0.740 






0.480 


0.580 


0.660 


Japan 


65 


0.397 ± 0.052 


0.631 






0.334 


0.631 


0.766 


(Japan — Non-scented) 


60 


0.397 ± 0.053 


0.623 






0.352 


0.639 


0.759 


Fragrant landraces scented 


35 


0.336 ± 0.055 


0.483 






0.000 


0.633 


0.738 


E13 SNP group 


29 


0.325 ± 0.055 


0.328 


0.000 


0.000 


0.000 


0.523 


0.709 


Splic-mutant group 


6 


0.217 ± 0.047 


0.444 


0.000 


0.000 


0.000 


0.000 


0.500 


Fragrant landraces non-scented 
















Non-scented group 


5 


0.253 ± 0.044 


0.000 






0.320 


0.480 


0.640 


Non-fragrant group 


25 


0.349 ± 0.048 


0.435 


0.480 


0.499 


0.570 


0.557 


0.509 



-: not genotyped. 
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Fig. 2. Phylogenetic tree based on NJ method. Bar showed genetic distance measured by 30 SSR genotypes. Nine clades (Clade 1 to 9) were 
shown as CI to C9. Bootstrap values were assessed with 1000 replicates and the scores were shown beside each clade. 



A. 



1296nt 



Up1136 (6/8bp-deletion) 



Badh2 



Upl098(+T) 



Up-201 (+G) 



5 0 0bp 

B. 



Splicing mutation s bp-deletion 

CCCATCGgtaccc Nipponbare AAAAAGATTATGGCTTGAGCTGC Nipponbare 
CCCATCGataccc Clade 8 AAAAA CTTCAGCTGC Clade 5 



-did D-D D-D-D-D D 

▲ 1 
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Fig. 3. Gene structure of BADH2 with upstream region. Sequence differences found in fragrant landraces and INDELs compared to reference 
sequence in Nipponbare and others were noted as solid triangles. The size of triangles corresponded to relative size of the INDELs. A. Three 
INDELs were recognized in upstream regions. B. Inside Badh 2, two functional SNP and deletions were recognized. Splicing junction at exonl- 
intronl, intron sequences were shown as small letters and g at splicing donor site was substituted into a showing bold letter. One known mutation 
was C to T substitution at exon 13. Three INDELs (solid triangles) and SSR (flag with solid circle) were detected. One of INDELs, 8bp deletion at 
exon 7 was carried by Koihonoka and some oversea landraces. 
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Table 3. Sequence substitution, insertion (+), and deletion (-) compared to Nipponbare sequence among representatives of different clades in 



Japanese fragrant landraces 









Upstream 








Tn tvr*~n 1 
111UU11 1 


Intron 2 


Tntvr^Ti A. 


tvi tivwi 1 
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0 
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-TT 
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40-nichi-wase 






























Iwaka 


Clade-3 


+AAATCTAG 


IT 
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T 
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G 


-TT 
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11 


T 


T 


G 


Kabashiko A 


Clade-7 


+AAATCTAG 
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G 


-TT 


c 


11 


T 


T 


G 


Nioi-mochi 2 


Clade-8 


+AAATCTAG 


+T 


A 


G 


T 


0 


A 


-TT 


c 


11 


C 


C 


G 


Hieri 


Clade-9 


+AAATCTAG 


+T 


A 


-G 


T 


0 


G 


-TT 


c 


11 


T 


T 


G 



change was detected only in Clade 8, which carried the 
same G to T (G/T) substitution in the splice donor site at the 
exon 1 -intron 1 junction in 'Nioi-mochi 2' (Fig. 4 and 
Table 3). The G/T substitution involves to a recognition site 
of Kpnl restriction enzyme. Non-fragrant rice such as 
'Nipponbare' and general fragrant rice can be digested with 
Kpn I, whereas fragrant landraces in Clade 8 can not be di- 
gested because of the substitution at the splicing donor site. 
Thus, a CAPS marker was developed to confirm the sequence 
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Jakou-ine 
Nipponbare 
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> 
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Fig. 4. Novel allele at Badh 2 locus. A. A single nucleotide substitu- 
tion was identified at exon-splicing junction, G to A, which is Kpn I 
resctriction enzyme recogniton site. Jakou-ine in Clade 8 carried the 
mutation type and reference genome Nipponbare, non-fragrant elite 
varieties carried a general type. B. PCR amplicon in Nioi-mochi 2 was 
not digested by enzyme digestion but others were. All fragrant land- 
races were applied and only six landraces in Clade 8 showed that 
their amplicons were not able to digest. 



substitution. All fragrant landraces were re-confirmed by 
Kpn I restriction digestion of amplicons including the rele- 
vant site (Fig. 4). As the SNP has not yet been reported, it 
was concluded that six varieties in Clade 8 carry the novel 
allele at the Badh2 locus. This has now been registered at 
DDBJ(AB821274). 

Since the six landraces belonging to Clade 8 shared the 
splicing donor site mutation, their morphological and physi- 
ological characteristics were compared. Three of them were 
non-glutinous and the remainder were glutinous (Fig. 5). 
Grain shape, presence or absence of awns, and colors in 
parts of the spikelets were polymorphic. 

INDEL polymorphism 

Based on the sequence data, six INDELs were detected 
among representatives of fragrant landraces (Fig. 3). These 



1 Jakou-ine 2 

M • t 

ii 1 Nioi-mochi 2 

Jakou-mochi w 

•I •# 



Fig. 5. Spikelets of six landraces in Clade 8. Nioi-mochi 1, Nioi- 
mochi 2, and Jakou-mochi represented glutinous endosperm and others 
were non-glutinous. Jakou-ine 1, Jakou-ine 2, Nioi-mochi 2, and 
Kamari were awned spiekelets. Nioi-mochi 1, Nioi-mochi 2, Kamari, 
and Jakou-mochi carried C and A alleles at C locus and A locus, re- 
spectively. Thus, their awns or apiculus part (top of the spikelet) and 
sterile glumes were colored. 
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Table 4. Haplotypes evaluated with inserts (+) and deletions (-) compared to 'Nipponbare' sequence and the frequencies of landraces and varie- 
ties among overseas fragrant carrying 8bp deletion or no-deletionat exon 8, E13 SNP, splicing donor site mutant (Splic-mut), and non-fragrant 
groups 



DNA variation in each haplotype No. of landraces/varieties 

Haplotype Upstream Inside BADH2 Overseas Japanese groups 





Upll36 


Up 1098 


5'UTR 


I2-986nt 


8 bp del 


non-8 bp 


E13 SNP 


Splic-mut 


Non-fragrant 


Indel-Hap 1 * 








-2 bp 


0 


0 


0 


0 


8 


Indel-Hap2 




+T 






0 


0 


0 


0 


7 


Indel-Hap3 




+T 




-2 bp 


0 


0 


0 


0 


3 


Indel-Hap4 


-6 bp 






-2 bp 


0 


3 


0 


0 


0 


Indel-Hap5 


-8 bp 


+T 


+3bp 


-2 bp 


1 


0 


0 


0 


0 


Indel-Hap6 


-8 bp 


+T 






2 


0 


3 


0 


7 


Indel-Hap7 


-8 bp 


+T 




-2 bp 


7 


4 


26 


6 


0 


Total 










10 


7 


29 


6 


25 



* 'Nipponbare' was included in Hap 1 . 



INDELs were detected through sequencing of the entire 
coding region including the 5 '-upstream region with Japa- 
nese fragrant landraces shown in Table 3, 'Delia', and 
'Meragome'. Only Delia showed the 3 -bp deletion at the 
5'-UTR, and subsequent screening of the deletion also indi- 
cated that it is unique in 'Delia'. The 8 -bp deletion at exon 7 
is known to cause loss of function and is found only in over- 
seas landraces in the C4-1 Clade and 'Delia' in the C4-2 
Clade (Table 4). 

In total, seven haplotypes were detected among overseas 
and Japanese groups (Table 4). Three INDELs, Up201, 
Upl098 and Upll36, were detected in the upstream region. 
Upll36 was a 6-bp or 8-bp deletion without reference to 
'Nipponbare'. However, seven non-fragrant varieties car- 
ried the 8-bp deletion for Upll36. Two single-nucleotide 
insertions were found in the Japanese fragrant group. One 
(Up201) could not be converted to an INDEL marker, but 
another (Up 1098) could, and this INDEL was also found in 
the non-fragrant group. 

There were two INDELs in transcribed regions, a 3 -bp 
deletion in the 5'-UTR and a 2-bp deletion within intron 2 
(1-986). The former deletion was found only in non- 
Japanese fragrant groups while the latter (1-986) was detect- 
ed in both non-fragrant and fragrant groups. These INDELs 
did not lead to loss of gene function. 

Allelism test 

'Koihonoka' is a recently released cultivar from Aomori 
prefecture, Japan, and carries the 8-bp deletion of exon 7 in 
Badh2 as a single causal mutation of fragrance (Fig. 6 and 
Table 5). This was confirmed using F 7 recombinant inbred 
lines (RILs) from the cross between the non-fragrant 
modern variety, 'Masshigura', and the fragrant variety, 
'Koihonoka'. One hundred sixty-five individuals were phe- 
notyped by sensory testing to confirm fragrance, and they 
segregated as 85 non-fragrant : 80 fragrant. This was com- 
patible with a 1 : 1 ratio (chi-squared test = 0.15), suggest- 
ing that the fragrance was controlled by a single gene. All 
progeny individuals were genotyped using a INDEL marker 
to show the 8-bp deletion. Complete linkage was obtained 



+ + - - + - + - + 




Fig. 6. A 8 bp deletion of exon 7 in BADH2 corresponding to pheno- 
types of fragrance in recombinant inbred lines of the F 7 progeny of 
'Masshigura' x 'Koihonoka'. Deletion type in 'Koihonokara' was in- 
herited from Basmati370. Phenotypes were segregated as a single 
gene. + and - meaning fragrance and non-fragrance by using leaf test 
with 2% KOH treatment. 

between fragrance and genotypes of Badh2. 

'Masshigura' was crossed with 'Koihonoka' and 'Jakou- 
ine 4' (Clade 8). In addition, other non-fragrant cultivars 
were crossed to 'Iwaka' (Clade 3), 'Kaori-ine' (Clade 3), 
and 'Jakou-mochi' (Clade 8) in order to test Mendelian seg- 
regation as a single gene. These F 2 populations segregated 
into non-fragrant and fragrant individuals in a 3 : 1 ratio. 
Thus, the fragrant trait was controlled by a single gene. 

Fj individuals of 'Koihonoka' x 'Iwaka', 'Koihonoka' x 
'Kaori-ine', 'Koihonoka' x 'Jakou-mochi', 'Iwaka' x 'Jakou- 
mochi', and 'Kaori-ine' * 'Jakou-mochi' were subjected to 
an allelism test. These Fj individuals were genotyped using 
SSR markers to confirm whether they were heterozygous 
as F^ Sensory testing suggested that all the F[ individuals 
and the self-pollinated progenies were fragrant. It was con- 
cluded that their fragrant alleles were controlled by the 
same single gene, Badh2, as recessive multiple alleles. 
There were at least two groups as Japanese landraces, 1) 
the E13 SNP group with the SNP in exon 13, as described 
above, and 2) the splicing mutant (Splic-mut) group carry- 
ing the G/T substitution in the splicing donor site at the 
exon 1 -intron 1 junction without the SNP at exon 13. Addi- 
tionally, there is another group included 'Koihonoka' and 
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Table 5. Allelism test of fragrant mutations among different landraces and 'Koihonoka' (modern variety) carrying the 8bp deletion in Badh2 in- 
heritted from Basmati 370 



F2/F7 populations 




No. of plants 


No. of plants scored for the 
fragrance 


Chi-test against a single 

'iPtJrpoatinTi f3 ■ 1 or 1 


gene 








Non-fragrance 


Fragrance 


K, ndTiiilntimic 
[^nirintint'a 


X C^laHp Iwaka 


10 


0 


10 


( All fpcpcsivp h/npl 

^rtii icccaMvc Lypc^ 




Koihonoka 


x Clade 3: Kaori-ine 


10 


0 


10 


(All recessive type) 




Koihonoka 


x Clade 8: Jakou-mochi 


10 


0 


10 


(All recessive type) 




Clade 3 : Iwaka 


x Clade 8: Jakou-mochi 


10 


0 


10 


(All recessive type) 




Clade 3 : Kaori-ine 


x Clade 8: Jakou-mochi 


10 


0 


10 


(All recessive type) 




Koihonoka 


x Masshigura 


50 


39 : 


11 


NS (0.5 < p) 




Clade 4: Narukosan-koutou 


x Clade 3: Iwaka 


50 


39 : 


11 


NS(0.5<p) 




Clade 3 : Kaori-ine 


x Clade 1 : Nioi-wase 1 


IS 


42 : 


8 


NS (O.Kp) 




Masshigura 


x Clade 8: Jakou-ine 4 


50 


38 : 


12 


NS (0.5 < p) (test for 3 : 


1 ratio) 


F 7 RIL 














Masshigura 


x Koihonoka 


165 


85 : 


80 


NS(0.5<p) (test fori : 


1 ratio) 



'Meragome' inherited the 8 -bp deletion in Badh2 from over- 
seas fragrant landraces. 

Genetic diversity around the Badh2 locus 

Heterozygosity was measured as the average score for 30 
SSR loci distributed over the whole genome (Fig. 7 and 
Supplemental Table 1). An additional 16 SSR markers were 
assayed within a 543 -kb stretch of Badh2 in order to identi- 
fy selective sweeps around the Badh2 locus (Supplemental 
Table 2). Fourteen of the additional SSRs were monomor- 
phic or nearly monomorphic as only one cultivar showed an 
allele different from the others. Such rare alleles are proba- 
bly due to a recent spontaneous mutation. Two additional 
loci, 45kup and 12kup, were relatively highly heterozygous, 
He ranging from 0.480 to 0.499, respectively, in the non- 
fragrant group. He scores at four loci, RM8264, SSRat7, 
RM1109, and RM3459, were compared with the data ob- 
tained at the two polymorphic loci. The splicing mutant 
group was monomorphic within at least a 151 -kb region of 
Badh2 and the El 3 SNP group was monomorphic within at 
least a 45-kb region. Splicing mutants were cultivated at the 
restricted area carried more wider selective sweep region. It 
may result from the recent origin of the novel mutation. 




Fig. 7. Selective sweep found around Badh 2 locus. Vertical line 
showed Heterozygosity measured by 30 SSR genotypes. Horizontal 
axis corresponded to each SSR locus. 



Discussion 

Fragrant rice is an economically important crop in terms of 
global trade. The FAO All Rice Price Index (2002-04 = 
100) indicates that among all types, fragrant rice varieties 
were priced more highly than others, as of January 2013. 
The causal gene for fragrance has been reported as Badh2 
(Kovachi et al. 2009). Ten multiple alleles have been found. 
One of them, Badh2.10, is shared between Taiwanese 
fragrant landraces and one Japanese landrace (Shiratama). 
As we have various Japanese fragrant landraces except 
for Shiratama, we characterized them at the molecular level 
in the present study to screen for novel mutations. When 
data of all SSR markers were analyzed together, we recog- 
nized nine different clades among modern cultivars and 
landraces, including non-fragrant landrace/varieties. The 
Japanese fragrant landraces were relatively highly diverse 
and formed Clades 1, 2, 3, 7, 8 and 9, which reflected their 
regional divergence. However, based on the Badh2 haplo- 
types, fragrant landraces comprised only two groups, one 
corresponding to Clade 8 and the other to Clades 2, 3, 7, and 
9. The latter group carried the SNP in exon 13 (E13 SNP 
group), which was reported by Kovach et al. (2009) as 
badh2.10. Although te El 3 SNP group shared the single al- 
lele, they showed higher structure which would result from 
adaptation to each environmental condition through out- 
crossing. The original lineage carying the SNP was sponta- 
neous ly outcrossed with other landraces for exchange of 
the genetic background. Although landraces in Clade 8 were 
fragrant, they did not carry the exon 13 SNP. Like the 
low-amylose-type Wx h allele at the Wx locus (Sano 1984) 
they shared a single mutation at the splicing donor site, 
which was the only mutation within the gene and upstream 
of it. INDEL haplotyping suggested that there were no 
INDELs related to the function of Badh2. Our allelism 
test also suggested that the splicing mutation was allelic 
to the 8-bp deletion allele shared with 'Basmati 370' and 
'Koihonoka', suggesting that the splicing donor site muta- 
tion results in fragrance. The exon 13 SNP mutation was 
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the most common among Japanese fragrant landraces while 
the remainder ( 1 7%) were found to carry the splicing donor 
site mutation whih has never been reported. 

Whole-genome genetic diversity estimated in terms of 
SSRs showed that the two groups of Japaneses fragrant land- 
races were not significantly different from the non-fragrant 
group, probably due to genetic exchange with the non- 
fragrant group. In contrast, estimation of genetic diversity 
with SSR markers flanking Badh2 showed that the non- 
fragrant group had higher genetic diversity in this region 
relative to the two groups of fragrant landraces, although the 
splicing mutant group carried a longer stretch of sequences 
without any polymorphisms than the E13 SNP group. Out- 
side this region, genetic diversity increased dramatically. 
The high diversity in fragrant landraces were also visible as 
morphological and physiological traits. Although the splic- 
ing donor site mutants shared a single mutation that had 
originated recently, they showed differences in endosperm 
traits, awns, and organ coloration, which had probably 
arisen due to outcrossing with different landraces. The lack 
of polymorphism around the Badh2 locus can be explained 
by the presence of a selective sweep similar to the situation 
in glutinous rice (Olsen and Purugannan 2002) where a re- 
duction of nucleotide variation spans a stretch of -250 kb 
centered on the Waxy locus (Olsen et al. 2006). 

Except for one Japanese fragrant landrace, Jakou-ine, 
variety names of the novel mutants were not listed in histor- 
ical records of Japanese fragrant landraces of the 17th to 
19th centuries, but do appear in records published in 1970 
(Arashi 1975). Anumber of landraces carrying the E13 SNP 
have been named independently, which have a single origin 
for the fragrant gene known as badh2.10. It was also culti- 
vated at Taiwan. However, the novel splicing mutation has 
not been found outside Japan. The distribution of the mu- 
tants is partly overlapped with that of fragrant soybean land- 
races in the same southern Tohoku area, one of which is 
known as the fragrant soybean, 'Dadacha-mame', a well 
known product that was recently found to carry a mutation 
in the Badh2 homolog, GmAMADHl (Arikit et al. 2011). 
Fukushima prefecture, which neighbors Yamagata, also has 
a fragrant soybean landrace. The origin of the splicing mu- 
tation in rice plants may be a recent event happened around 
the area. No variety sharing the splicing mutation was found 
in other areas. Preference of local people allowed the selec- 
tion as an independent lineage to create novel fragarnt 
crops. Probable conscious selection such as this has been 
noted for Indian corn, which contrasts with more recent 
negative selection for traits such as red pericarp in Japanese 
rice (Arashi 1976, Ishikawa et al. 2002a). Because the local- 
ity, landraces in the splicing mutant group was only com- 
posed from Tm-japonica, but not from Tx-japonica due to 
the local origin and genetic background cultivated at the 
area (Supplemental Table 5). 

The quantity of 2AP is influenced by environmental con- 
ditions and may be affected by seasonal and storage factors 
(Itani 2004). Multiple alleles may also alter the activity of 



enzymes that pass precursors down the biosynthetic path- 
way. The amino acid change originating from the E13 SNP 
group does not lead to complete loss of function in the en- 
zyme BADH2, and probably explains the lower concentra- 
tion of 2AP (Kovach et al. 2009) in this group or fragrant 
landraces relative to other groups or fragrant rice. Variation 
in the level of fragrance is a useful breeding target, as fra- 
grant varieties differing in this trait can be offered to con- 
sumers. Furthermore, these landraces have acquired various 
genetic constituents through outcrossing. Fragrant landraces 
without any aroma would be segregating descendants of 
such spontaenous crossing. In addition, higher structure of 
fragrant landraces sharing the single gene, badh2.10 in El 3 
SNP group, resulted from the spontaenous outcrossing. This 
lead them to adapt to various environmental condition. Such 
spontaenous outcrossing, however, has never been allowed 
in modern varieties. Because original seeds of modern vari- 
eties are strictly controlled to exclude outcrossed seeds from 
nursery fields in which seeds are propagated for farmers. 
Captured genetic constituents from different landraces 
might be lost from current genetic resources. In this sence, 
fragrant landraces are a genetic reservoir that could allow 
rice breeders to find unexpected traits such as tolerance to 
biotic and abiotic stress, which had been lost in elite varie- 
ties. Some presumed tolerance could be recognized under 
varietal names (Arashi 1975). Parts of that would be benefi- 
cial for modern varieties to acquire further valuable traits. 
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